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Self-propelled colloidal objects, such as motile bacteria or synthetic microswimmers, have mi¬ 
croscopically irreversible individual dynamics - a feature they share with all living systems. The 
incoherent behavior of individual swimmers can then be harnessed (or “rectified”) by microfluidic 
devices that create systematic motions impossible in equilibrium. Here we present a computational 
proof-of-concept study, showing that such active rectification devices might be created directly from 
an unstructured “primordial soup” of light-controlled motile particles, solely by using spatially mod¬ 
ulated illumination to control their local propulsion speed. Alongside both microscopic irreversibility 
and speed modulation, our mechanism requires spatial symmetry breaking, such as a chevron light 
pattern, and strong interactions between particles, such as volume exclusion causing a collisional 
slow-down at high density. Taken together, we show how these four factors create a novel, many- 
body rectification mechanism. Our work suggests that standard spatial-light-modulator technology 
might allow the programmable, light-induced self-assembly of active rectification devices from an 
unstructured particle bath. 


One-sentence summary: Self-propelled particles that swim in response to light can self-assemble 
microfluidic rectification devices under non-uniform illumination. 


INTRODUCTION 

The last decade has seen a surge of interest in so-called 
“active” materials, whose constituent building blocks vi¬ 
olate time-reversal symmetry on the microscopic scale 
by continuously converting fuel into motion . These 
building blocks can be natural, such as swimming bacte¬ 
ria or algae [3], or synthetic, such as self-phoretic Janus 
swimmers Apart from deepening our understand¬ 

ing of nonequilibrium systems, and in particular about 
how and when concepts from equilibrium thermodynam¬ 
ics can (and cannot) be applied to such systems [6j[7], a 
strong motivation behind active matter research is that 
activity can be harnessed to create devices and induce 
phenomena that would be impossible in thermodynamic 
equilibrium. Examples include flow of rotor particles 
round a circuit [8j, steady rotation of a gear wheel in a 
bacterial bath 0 HD], and pumping of bacteria between 
chambers by “funnel gates” m- 

In this study, we show how irreversible motion in a 
particular class of active matter - light-controlled self- 
propelled particles - can be used for self-assembling ac¬ 
tive “rectification devices” and circuits directly from an 
unstructured active-particle suspension using only spa¬ 
tially controlled (but possibly incoherent) lighting. Such 
light-controlled motility has been experimentally demon¬ 
strated several times, both in synthetic systems HMDS 
and natural ones (strains of Escherichia coli bacteria 
whose motility requires illumination) [16]. The possibil¬ 
ity of such light-controlled self-assembly is far from ob¬ 
vious, since light intensity is, for our purposes, a scalar 
rather than a vector field: it controls the speed of par¬ 


ticles, but not their orientation; therefore, uniform illu¬ 
mination does not itself create directional motion. This 
differs fundamentally from magnetotaxis, for example, 
in which a uniform vector field directs motion by align¬ 
ing the propulsion direction of particles PH- Because 
of the ubiquity of spatial-light-modulator (SLM) tech¬ 
nology [18], the ability to create active rectification de¬ 
vices using light-induced self-assembly could have many 
advantages over conventional microfluidic approaches to 
rectification mmmm- In particular, it would allow 
real-time device reprogramming rather than mechanical 
interchange of fixed microfluidic elements. Our proposed 
strategy hinges on the tendency of motile particles to 
accumulate in regions where they move slowly [221124] . 
This allows one to use light to sculpt regions of high and 
low particle density. So long as the particles interact re¬ 
pulsively, the dense regions might be used to create, in 
effect, microfluidic obstacles of a type previously used for 
rectification m • We first develop this appealingly sim¬ 
ple avenue before turning to other approaches that prove 
more efficient. 

The slowing-induced accumulation of motile particles 
is much like the crowding of pedestrians who slow down 
in front of a shop window on a busy street [6]. Within a 
wide class of models 24] . it is quantified as p( r) oc l/v(r ), 
where p{ r) is the steady-state mean particle density at 
position r and v(r) is the local swim speed 22]. This 
result has no counterpart for nonmotile (e.g., Brownian) 
particles in thermal equilibrium, whose steady-state den¬ 
sity is governed solely by the Boltzmann distribution, and 
whose speed statistics remain completely independent of 
their position. It applies rigorously for noninteracting 
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particles [22[ [24] and generalizes as p( r) oc l/v(r , p) to 
cases where interactions cause the average local swim 
speed v(r,p) to differ from the local value for an iso¬ 
lated particle, v (r), which we assume to be illumination- 
controlled [23j '24b 


RESULTS AND DISCUSSION 

Light-controlled accumulation of active particles 

For our proof-of-concept simulations we deploy a well- 
studied model of “active Brownian particles” in two spa¬ 
tial dimensions. Each particle swims with a constant 
speed v along a particle-fixed axis that rotates by au¬ 
tonomous angular diffusion; this motion is influenced 
by interparticle collisions (see Methods). It is known 
that at fixed, spatially uniform v , collisions cause a lin¬ 
ear dependence of the average swim speed on density: 
F(0) = ^(l — 0/0*) with 0 = pira 2 /4 the volume fraction 
of particles in two dimensions, a the particle diameter, 
and 0* « 0.93 a near-close-packed volume fraction that 
depends on interaction details [23 26]. Hence, in this sys¬ 
tem one might expect that if a spatial speed-modulation 
pattern v(r) is imposed, 0(r) should everywhere obey 
0 oc l/(v(l — 0/0*)). This prediction is tested in Fig. [l] 
for a modulation pattern comprising a lattice of 16 cir¬ 
cular traps, each with a different uniform nominal speed 
Vi = vq(1 — i/16), i = 1 ,..., 16, in its interior and a larger 
speed vo between the traps. (We use periodic boundary 
conditions.) The agreement is excellent, except at the 
highest densities where, with the chosen interactions (see 
Methods), particles show significant overlap, invalidating 
the assumed linear dependence of v on 0 [25 . 


Self-rectification by funnel gates 

The above results demonstrate that for motile particles 
with repulsions, just as for noninteracting ones, the den¬ 
sity pattern p(r) can be rather precisely manipulated by 
spatial modulation of the one-particle swim speed i?(r). 
But Fig.[l]does not demonstrate rectification: that would 
require either a continuous steady-state particle current, 
or, if the boundary conditions prevent such a current, two 
areas of the same v to acquire different densities. Only 
then can one claim to have converted incoherent particle 
dynamics, which by itself is responsible for the p oc 1/v 
behavior, into a systematic pumping effect. 

A simple strategy for achieving such rectification is 
to create a light-guided density pattern that mimics mi¬ 
crofluidic devices. Such devices generally operate at the 
one-particle level: a series of obstacles is created so that 
each particle is geometrically biased toward a nonzero 
average motion. However, when the obstacles are made 
of other particles this becomes an inherently many-body 



v(l-4>/<f>*)/vo, v/v 0 


Figure 1. Time-averaged local volume fraction 0(r) of in¬ 
teracting and noninteracting particles in a system with spa¬ 
tially patterned one-particle swim speed u(r), where v — Vi — 
uo(l — i/16) in the ith well and v = Vo between the wells. The 
solid lines show our predictions with a prefactor determined 
by least-squares fits, where we neglected the three highest 
density (0 > 0.5) data points for the interacting particles. 
Inset: Spatially resolved plot of 0(r) for the interacting par¬ 
ticles; note the logarithmic scale on the color bar. 


mechanism that computer simulations can illuminate. In 
an elegant experimental demonstration of one-particle 
rectification, Galajda et al. El established that a wall 
of microfabricated chevron-shaped obstacles or “funnels” 
will rectify a suspension of swimming Escherichia coli 
bacteria. Asymmetric pumping across such a wall be¬ 
tween two chambers of the same v can yield more than 
a factor-two density difference. The effect of additional 
layers of walls is multiplicative on the density ratio m- 

Inspired by these experiments, we now simulate a set 
of chevron-shaped “dark” regions, in which the nomi¬ 
nal swim speed u(r) is set to zero. As can be seen in 
Fig. 1 and Movie SI, this pattern rapidly guides the 
self-assembly of a set of funnel-like obstructions which 
in turn cause rectification. The mechanism is almost 
the same as in Refs. mm- the trajectories of motile 
particles will tend to trace the funnel wall, making the 
probabilities of passing through the funnel from either 
side unequal. However, because our funnels assemble re¬ 
versibly from the motile particles themselves, such rec¬ 
tification requires not only time-irreversible trajectories 
(wall-tracing) and broken reflection symmetry of the 
modulation pattern PS EH EE] but also repulsive inter¬ 
particle interactions (Fig. |2|. Unlike conventional micro- 
fabricated structures, ours could in principle be disman¬ 
tled and re-assembled into a completely different struc¬ 
ture, simply by reprogramming the light intensity pat¬ 
tern to give a different f(r). Crucially, the mechanism 
requires no external or interparticle forces that align the 
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Figure 2. (a) Time-averaged particle packing fraction </>(r) at 
steady state in a system of self-assembled funnels. Note the 
logarithmic scale on the color bar. (b) Magnification of one of 
the funnels, showing the self-assembly process at three differ¬ 
ent times starting from a homogeneous particle distribution; 
tlj is the Lennard-Jones time (see Methods), (c) Black curve: 
The packing fraction 0 from (a), but now as a function of x 
and averaged over y\ the pumping effect (rectification) is ob¬ 
vious. Red curve: For noninteracting particles, rectification 
is not possible. 


ria preferentially transport passive tracer colloids from 
the shallow side and deposit them on the steep side of 
the barrier, leading to a significant steady-state accu¬ 
mulation. Loosely inspired by this, we abandoned the 
chevron design, and imposed a simpler, one-dimensional 
sawtooth-shaped speed profile v(x) shown in Fig. [ 3 ] This 
leads to the self-assembly of an asymmetric barrier with 
a distinctive density profile p(x). 



Figure 3. Steady-state packing fraction (f){x) (solid curves) 
in a system with a modulated speed v(x) composed of two 
opposing saw-tooth profiles (dashed curves), obtained from 
particle-resolved simulations and our continuum model. For 
numerical reasons, the discontinuities of v{x) are smeared out 
over 10cr in the continuum calculations. The inset shows (j>(x) 
with a magnified y scale, clearly demonstrating the rectifica¬ 
tion. 


swimming directions of the particles. 

As a first proof of principle, these are striking results. 
However, the efficiency of our device (measured as the ra¬ 
tio between the densities in the inner and the outer com¬ 
partments) is modest: the density ratio induced by three 
pairs of funnel walls is about 1.6 (see Fig. [ 2 J 3 ) , compared 
to about 2.8 seen in Escherichia coli experiments using 
a single microfabricated funnel wall El and in simula¬ 
tions of rod-shaped particles in a similar geometry m- 
As shown in the Supplementary Material, this is largely 
because our self-assembled many-particle funnels cannot 
prevent particles from “leaking” backwards against the 
pumping flux, limiting their pumping efficiency. 

Self-rectification by one-dimensional speed profiles 

A contrasting single-particle rectification strategy de¬ 
ploys microfabricated wedge-shaped solid barriers with 
one steep and one shallow side. Such barriers have been 
shown experimentally [201 ED to make swimming bacte- 


Remarkably, this extremely simple device shows a den¬ 
sity ratio of almost 2 and thus a strong rectification (see 
inset of Fig. [3|. Note that particles must pass through, 
but not over or around, the barrier so that the “leakage” 
encountered earlier is being exploited here. To study 
this rectification, we have extended a continuum model, 
involving local density and orientation fields coupled in 
a standard manner [29], to the case where the swim 
speed depends on both position and density (see Meth¬ 
ods). The rectification is found to arise from a polariza¬ 
tion (i.e., alignment of particle orientations) on the dense 
(steep) side of the barrier; here the average swimming di¬ 
rection of the particles is oriented into the barrier (see the 
Supplementary Material for details). This polarization is 
a straightforward consequence of the irreversible particle 
motion, which makes particles that swim towards a solid 
object such as a wall become “trapped” there until their 
swimming direction has rotationally diffused enough to 
let them escape, creating an average particle orientation 
pointing into the object. The polarization in turn creates 
a mechanical pressure that drives a particle flux from the 
steep side of the barrier towards the shallow side. In more 
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mathematical terms, rectification rests on the presence of 
a density-gradient-driven pumping effect that is nonzero 
only in a small region close to the barrier. Flux bal¬ 
ance then requires a change in the uniform steady-state 
particle densities in the bulk regions far to the left and 
to the right of the barrier, where density gradients, and 
thus the pumping effect, approach zero. In contrast to 
equilibrium systems, here local gradient effects cause rec¬ 
tification across an interface between bulk regions, rather 
than merely shifting the interfacial tension [ 501 . 


Self-assembled active circuits 

By using aligned (instead of opposing) barriers and ap¬ 
plying suitable boundary conditions, either of the above 
devices can be reconfigured to create an “active circuit”, 
where motile particles flow endlessly around a closed loop 
at steady state. In Ref. [8j macroscopic circulation was 
achieved with field-driven autorotors with hydrodynamic 
alignment interactions. In Fig. [4] we give computational 
proof-of-principle for a simpler design that employs the 
speed modulation from Fig. [3j but now with two aligned 
saw teeth. In this example, colloids are confined to an 
annulus by a surrounding dark region; thus, both the 
“circuit” itself and the rectifying elements arise solely by 
light-guided self-assembly. 



Simulation Theory 


Figure 4. (a) Particle-resolved simulations of an annular cir¬ 
cuit self-assembled through a prescribed speed profile u(r) 
with two similar saw teeth along the annulus and a very 
low speed outside of the annulus. The resulting dense par¬ 
ticle packing outside of the annulus (not shown for clarity) 
confines the motile particles to the circuit. The color map 
shows the time-averaged steady-state local packing fraction 
0(r) and vectors denote the rescaled local particle current 
J(r), with their colors indicating the absolute current density 
J(r) =s | J (r) |. The logarithmic color bar applies both to the 
density plot and to the vector field plot, (b) Corresponding 
plot of the rescaled particle density p(r) oc 0(r) and particle 
current J(r) obtained from our continuum model. 


CONCLUSION 

In conclusion, we have given proof-of-principle evi¬ 
dence, using computer simulations and theoretical con¬ 
siderations, that self-assembled active rectification de¬ 
vices for motile particles can be created solely through 
spatial modulation of the local propulsion speed f(r). 
This space-dependent propulsion speed can in turn be re¬ 
alized by a space-dependent illumination HMMMS] 
and hence rendered programmable through use of SLMs. 
A key element of the light-guided assembly strategy pre¬ 
sented here is that the illumination supplies the parti¬ 
cles only with scalar information, namely the instruction 
to swim at a certain speed. In some experimental sys¬ 
tems the light elegantly also supplies the energy source 
for swimming [13], but this is not always the case fl2j , 
and is not an essential ingredient in our proposed strat¬ 
egy. Crucially, in our mechanism the optical field does 
not exert directional forces directly on the particles. This 
differs fundamentally from, and requires much lower light 
intensities than, the manipulation of particles by optical 
tweezers; as an example, the light-propelled particles of 
Ref. [13] require light intensities of < 5 /iW//im 2 , com¬ 
pared to intensities of 10 2 — 10 4 /iW//im 2 typically used 
in optical tweezers dU El]. Our mechanism can be eas¬ 
ily applied to large-scale systems containing many par¬ 
ticles, circumventing several disadvantages of traditional 
micro-manipulation techniques m- such as significant 
heating-up of the system and the need for complicated 
optical equipment. Indeed, although use of SLMs would 
certainly be appealing, there is no requirement that the 
light field even be coherent. 

Our work points towards a new technological opportu¬ 
nity in which light-guided self-assembly is not only used 
to make crystals and other aggregates [T2l fT3l [T5] but 
also to create programmable, active rectification devices 
from an initially structureless “primordial soup” of self- 
propelled particles. The capacity to do this stems from 
an interplay between microscopically irreversible parti¬ 
cle motion, spatially asymmetric light modulation, and 
interparticle repulsions. These findings may also be rele¬ 
vant for controlling mass transport in biological systems, 
steering pattern formation in future applications of active 
particles in materials science, and the wider deployment 
of active particles in colloidal self-assembly. 

METHODS 
Particle-based model 

Our particle-resolved computational studies use a sim¬ 
ple model of swimming microorganisms (like motile bac¬ 
teria and protozoa) and artificial microswimmers, known 
as “active Brownian particles”. Such particles exhibit 
persistent propulsion in a particle-fixed direction which 
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rotates through thermal rotational diffusion. This model 
has been shown to capture many essential features of ac¬ 
tive matter systems such as bacterial suspensions, while 
still allowing for both the derivation of exact results from 
first principles and computational studies of systems con¬ 
taining millions of particles. The active colloids are mod¬ 
eled as spherical particles with diameter a moving on 
a two-dimensional substrate. They interact repulsively 
through a truncated and shifted Lennard-Jones potential 
U(r) = 4 s{{(j/r) 12 — (a/r) 6 )-\-s with the center-to-center 
particle distance r and an upper cut-off at r = 2 1 / 6 cr, be¬ 
yond which U = 0. e determines the interaction strength 
and /3 = 1/(£;bT) is the inverse thermal energy. The 
motion of these particles is described by the overdamped 
Langevin equations 

Ti = u(rj)u(<9j) + /3D T Fi({rj}) + ^2EP A T , (1) 

Qi = a/2Z)rAr , (2) 

where r i(t) is the position and 0i(t) is the orientation of 
the ith particle at time t. In our particle-resolved Brown¬ 
ian dynamics simulations we solved these Langevin equa¬ 
tions numerically using the LAMMPS molecular dynam¬ 
ics package [331. Here, v(r) is the nominal swim speed of 
an active particle at position r, u (6) = (cos(0), sin(#)) T is 
an orientational unit vector that denotes the propulsion 
direction of a particle with orientation #, and Fi({rj}) 
is the total conservative force on particle i, which re¬ 
sults from the interactions with other particles and thus 
depends on U (r). Dt and Dr = 3Z>r/cr 2 denote the con¬ 
stant translational and rotational diffusion coefficients of 
the particles, respectively, and Ax(t) and Ar(£) are zero- 
mean, unit-variance Gaussian white noise terms. The 
natural time unit in this model is the Lennard-Jones time 
tlj = a 2 /(s/3Dt)- Further technical details about the 
particle-resolved simulations can be found in the Supple¬ 
mentary Material. 

Continuum model 

In our continuum simulations we solved the following 
dynamic equations for the local particle density p(r, t) 
and the local polarization P(r, t) numerically: 

P = -V-(vW)+D t V 2 p, (3) 

W = - ^V(vp) - 7 iW - 7 2 W 2 W + A;V 2 W (4) 

- u>i(W-V)W + w 2 S7(W 2 ) . 

Here, W = pP is a weighted local polarization with 
modulus W = |W| and 71 , 72 , uq, and uq are pa¬ 
rameters. This continuum model can be obtained by 
an appropriate adaptation of the derivation of a mean- 
field theory for a suspension of active particles in Ref. 
[29]. The average swim speed u(r, p) of an active parti¬ 
cle at position r in an active suspension of density p is 


known from particle-resolved simulations |25[ (26] to be 
linearly decreasing with p. Nevertheless, in our numeri¬ 
cal calculations based on the continuum model we used 
the expression u(r, p) = u(r) exp(— <ap(r, £)), where the 
parameter a determines how quickly the average local 
swim speed decreases with density. This choice increases 
numerical stability without changing the qualitative con¬ 
clusions drawn here. 

The parameters 71 and 72 are positive and ensure that, 
in steady state and in absence of any density or veloc¬ 
ity gradients, W vanishes everywhere, k is an effec¬ 
tive diffusion coefficient for W and ensures continuity 
of this field. The term proportional to w\ encodes self- 
advection, whereas the term proportional to uq can be 
viewed as an “active pressure” term, as it can be written 
as —Vp with p = —W 2 W 2 . Both these terms generically 
arise when coarse graining models of self-propelled parti¬ 
cles j29ll34ff36] : the ensuing analysis suggests uq > 0 [29]. 
Furthermore, if one considers W as a fast variable and 
performs a quasi-stationary approximation [24], the W 2 
term leads to a term proportional to — rc 2 V 2 ((V(up)) 2 ) 
in Eq. its form is similar to that of the leading-order 
active contribution in the minimal field theory for active- 
particle phase separation (“Active Model B”) [30] . Fur¬ 
ther technical details about the numerical solution of the 
continuum model can be found in the Supplementary Ma¬ 
terial. 


SUPPLEMENTARY MATERIAL 

Supplementary material for this article is available at 
[URL to be inserted by the editor}. 

Movie SI: Movie showing the self-assembly process of 
a single chevron-shaped obstacle. 

Fig. SI: Steady-state density profiles for funnel-wall sys¬ 
tems with the funnel particles either mobile or fixed. 
Fig. S2: Local polarization around a sawtooth-shaped 
speed profile, obtained from particle-resolved simulations 
and the continuum model. 

Fig. S3: Polarization profiles obtained from the contin¬ 
uum model with W 2 = 0 and W 2 > 0 . 

Fig. S4: Propulsion speed profiles used in simulations 
of “active circuits”. 
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